Generation of vortex dipoles in superfluid Fermi gas in BCS limit 
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We theoretically investigate the generation of the vortex dipoles in superfluid Fermi gas in the 
BCS limit. The vortex dipoles are generated in superfluid either by moving an obstacle above a 
critical speed or due to the decay of the shock waves obtained on the sudden mixing of two superfluid 
fragments. We observe that in pancake-shaped traps, the shock waves can lead to the formation of 
density ripples, which decay into vortex dipoles due to the onset of snake instability. 



PACS numbers: 03.75.Ss, 03.75.Kk, 03.75.Lm 
I. INTRODUCTION 

One of the unique feature of the superfluids is that 
these can support only quantized vortices. In this 
context, the experimental observation of the vortex 
lattice across the Bardeen- Cooper- SchriefTer to Bose- 
Einstein condensate (BCS-BEC) cross-over has unequiv- 
ocally proved the superfluid nature of these systems pQ. 
In the experiment, vortices were created by rotating the 
trapped superfluid Fermi gas using a blue detuned laser 
beam. In a recent experiment [2 with pancake-shaped 
BEC, a pair of vortices with opposite circulation were 
generated by moving an obstacle potential above a criti- 
cal speed, which is equal to the fraction of the local sound 
speed. This vortex- ant i vortex pair is nothing but two 
dimensional analogue of vortex rings. In Bose-Einstein 
condensates (BECs), vortex rings have been experimen- 
tally generated from the dark solitons, which undergo 
long wavelength transverse instability called snake insta- 
bility [3 . Vortices and vortex dipoles can also be gen- 
erated by merging and interfering of multiple trapped 
BECs [4 . The decay of shock waves into vortex dipoles 
has also been observed in BECs [5]. The decay is medi- 
ated by the formation of solitons which decay via snake 
instability. In the non-linear Kerr-like medium, the prop- 
agation, non-linear response, and collisions between the 
dispersive superfluid-like shock waves have been experi- 
mentally observed [6.. Shock waves have also been gener- 
ated by merging or splitting BECs non-adiabatically [7]. 
The same method has been used to generate the shock 
waves in superfluid Fermi gas at unitarity [8]. Unlike 
BECs, the soliton mediated decay of shock waves into 
vortex dipoles has not been unambiguously detected in 
superfluid Fermi gases. In a recent work [9 , formation 
of shock waves and domain walls in strongly interact- 
ing Fermi superfluid has been investigated theoretically, 
and the possibility of snake instability was conjectured in 
sufficiently wide traps. The formation and dynamics of 
sub- and supersonic shock waves, using zero temperature 
equations of generalized superfluid hydrodynamics, have 
also been studied in unitary Fermi gas [TO] . 

Although vortex dipoles and vortex rings are yet to 
be experimentally realized in superfluid Fermi gases, it is 
well established that vortex antivortex pairs play very im- 



portant roles in superfluid turbulence [TT] and Berezinskii 
Kosterlitz Thouless (BKT) phase transition [12] [13]. Re- 
cently using local extension of the time-dependent den- 
sity functional theory, real time dynamics of quantized 
vortex rings in unitary Fermi gas in cylindrical traps was 
theoretically investigated in Ref. [14]. 

The advent of ultracold fermions in optical lattices [15] , 
which are almost pure realization of Hubbard model [16 , 
has lead to the flurry of research investigations. These 
systems are now routinely used, both by experimental- 
ists and theorists, to understand the strongly correlated 
systems. In this context, Mott insulator of fermionic 
atoms in an optical lattice has already been experimen- 
tally realized [17j[T8]. On the theoretical front, dynam- 
ical mean field theory (DMFT) [19] has been success- 
fully used to study these strongly correlated systems in 
three dimensions [20l [21]. The DMFT neglects the non- 
local correlations, which can no longer be done in lower 
dimensions and frustrated systems. For these systems, 
cluster extension of DMFT like dynamical cluster ap- 
proximation (DC A) and cellular dynamical mean- field 
theory (CDMFT) have been used to study these sys- 
tems [22H24] . A key recent development in the field, has 
been the experimental realization of spin-orbit coupling 
(SOC) [23 [26]. The interplay between SOC and inter- 
atomic interactions leads to many interesting phenomena 
[27] . Recently, finite temperature phase diagram of two- 
component atomic Fermi gas with population imbalance 
in the presence of Rashba spin orbit coupling has been 
studied in Ref. [28]. 

In the present work, we numerically study the gener- 
ation of vortex dipoles in an oblate superfluid Fermi gas 
(SFG) in BCS limit. We employ two methods, namely 

(a) vortex dipole generation by a moving obstacle and 

(b) soliton mediated decay of shock waves, to this end. 
At T = OK, when the normal component in the super- 
fluid tends to zero, a theoretical approach based on the 
Gallilei- invariant density functional theory has been de- 
veloped in Refs. [29H37] to study the superfluid Fermi 
gases. This theoretical approach allows one to write a 
non-linear Schrodinger equation (NLSE) for the BCS su- 
perfluid that leads to the same superfluid density as the 
original many body fermion system [30] . For the applica- 
bility of this NLSE for the BCS superfluid, the character- 
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istic wavelength of the phenomenon under study must be 
larger than the healing length. The NLSE has been used 
to study the collective excitations [30j [32], free expan- 
sion of the superfl uid Fermi gas [33] , superfluid-insulator 
transition [38] , and solitons [35], [36] . 

The paper is organized as follows. In Sec. [II] we study 
the evolution of the BCS superfluid with a repulsive ob- 
stacle moving across it. In Sec. [ITT] we study the gen- 
eration of the shock wave by non-adiabatically merging 
two superfluid fragments followed by the decay of shock 
waves into vortex dipoles. This is followed by conclusions 
in Sec. HVl 



here i/j = ^(r,£). In the present work, we consider 
pancake-shaped trapping potential 



V(r,t) 



V" {x 2 +y 2 + a 2 z 2 ) + y ohs ^ t ^ (3) 



where uo is the radial trapping frequency, a ^> 1 is the 
ratio of axial to radial trapping frequency, and V h s is the 
obstacle potential which can be attractive or repulsive. 
We consider SFG of 40 K with N = 10 3 , u = 10Hz, a = 
10, and 



V ohs (r,t) = Vbexp < -2 



[(x - x (t)) 2 + (y - y (t)) 2 



II. GENERATION OF THE VORTEX DIPOLES 
BY THE REPULSIVE POTENTIALS 



Barring normalization factors, the order parameter for 
the superfluid bosonic and fermionic systems are 



S B (r,t) = Mr,*)), 
S F (r,t) = (^(r,£)^ t (r,t)), 



(1) 



respectively, where ^(r, £) and ^(r, t) with a =f or ^ 
are the bosonic and fermionic annihilation field opera- 
tors, respectively. In case of bosons, the order parame- 
ter is the wavefunction of the macroscopically occupied 
single particle state. Similarly, the order parameter can 
be considered as the wavefunction of the macroscopically 
occupied two particle state for Fermi superfluids. The 
order parameter is normalized to total number of con- 
densed bosons [39] or condensed Cooper pairs [401143] . 
In the BCS limit, the s-wave scattering length between 
two hyperfine states of the Fermi gas a — > —0. In this 
limit, the superfluid Fermi gas, consisting of equal num- 
ber of two components, at T = K is described by mean 
field equation 



4m 



2h 2 



m Tj 



{^ 2 ) 2 ^(r,t)\^ + V(r,t) 
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= ih 



dt ' 



(2) 



where ip(r,t) is complex order parameter for the SFG, 
m is the atomic mass of the fermionic species, m p = 
2m is the mass of the fermionic pair, and V(r) is the 
trapping potential. In the aforementioned equation, 
ip(r,t) is normalized to the number of fermions iV, i.e., 
J ^(r, t)\ 2 dr = N; this normalization is different from 
the normalization of the Sp(r, t) [39 , 42 , 43 . In terms of 
^(r, t) the energy of the SFG is [29} [34] 
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here (xq (t) , yo (t) , 0) is the instantaneous location of the 
Gaussian obstacle potential with 1/e 2 width equal to wq. 
In Ref. [44], effective two dimensional (2D) equations for 
the pancake-shaped Fermi superfluids in both the BCS 
and unitary limit (a — >• — oo) have been proposed. These 
2D equations provide a good approximation to the real 
three dimensional (3D) systems. In the present work, in- 
stead of using these 2D equations, we consider the full 
three dimensional equation and solve it numerically us- 
ing the split time step Crank-Nicolson method [45 . We 
can rewrite the Eq. [2] in scaled units using the transfor- 
mations 



^(r,f) = 



t t'uj- 1 , 
N^r^t) 
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&OSC 
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where the primed quantities are in scaled units and a osc = 
^Jh/{m p u)) is the oscillator length. After dropping the 
primes, the scaled NLSE describing the superfluid Fermi 
gas in BCS limit is 



2(37r 2 AT) 2 / 3 |0| 4 / 3 + V(r, t) 



0,t) 



. d<P(r,t) 
1 dt ' 
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We consider the repulsive Gaussian obstacle potentials 
with Vo = 93/kj, xq = — 6a OS c 7 ^o(O) = 10.0/im and move 
it along x-axis. While moving the obstacle we contin- 
uously decrease the strength of obstacle potential such 
that it vanishes at xq = 6a osc . As is expected, we find 
that there is no generation of the vortices if the obsta- 
cle is moved below the critical speed. For example, Fig. 
[I] shows the evolution of the superfluid density when the 
obstacle is moved with the sub-critical speed of 150/im/s. 
On the other hand, when the obstacle is moved with 
a speed greater than the critical speed, vortex dipoles 
are produced. The generation of a vortex dipole for the 
obstacle potential moving with a super-critical speed of 
175/im/s is shown in Fig. [5J The critical speed for vortex 
nucleation is a fraction of the speed of the local speed of 
sound. In case homogeneous BCS superfluid, the speed 
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FIG. 1: (Color online) Evolution of the superfluid Fermi gas 
in BCS limit when the obstacle is moving with the sub-critical 
speed of 150 fim/s. Images with labels (a), (b), (c) are the 
superfluid densities |0(r,£)| 2 in units of a~ s 3 c on xy plane at 
t = ms, t — 143 ms, and t — 286 ms respectively. 



of the sound [46] c = v F /V%, where v F = 2(3tt 2 7V|^| 2 ) 1/3 
(in scaled units) is the Fermi velocity. For the parameters 
considered in the present work, the maximum speed of 
sound at the trap center is ~ 1465/im/s. This is only an 
approximate value since the presence of the trapping po- 
tential changes the speed of the sound from its value for 
homogeneous case [47] . When the obstacle moves across 
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FIG. 2: (Color online) Evolution of the superfluid Fermi gas in 
BCS limit when the obstacle is moving with the supercritical 
speed of 175 /xm/s. Images with labels (a), (b), and (c) show 
the superfluid densities |0(r, t)| 2 in units of a~^ c on xy plane 
at t = ms, t = 48 ms, and t — 127 ms respectively. 



the superfluid with the sub-critical speed, it can still gen- 
erate sound waves. The generation of the sound waves 
can be inferred from the variation of superfluid kinetic 
energy with time. The total kinetic energy KE of the 
superfluid in scaled units is 



KE 



dr 



dr. 



KE a + KE S 



(6) 



where 9 is the phase of superfluid wavefunction <\>. The 
first energy term KE q is the quantum pressure energy 
and the second KE S is the energy arising from the su- 
perfluid velocity. After the obstacle potential has be- 
come zero, the total energy almost becomes constant as 
is shown in Fig. |3|a). At the end of the time evolution, 
non-zero value of KE S is mainly due to the presence of 
swirls, i.e. due to the vortices, and sound propagation; 
Fig. [3^b) shows the change in KE S with time. In case 
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FIG. 3: (Color online) Variation of (a) E and (b) KE S with 
time. The dashed dark green curve shows the variation when 
the obstacle is moving with the speed of 175 /xm/s. The solid 
orange curve shows the variation when the obstacle is moving 
with the speed of 150 /xm/s. 

when the obstacle moves with the speed of 175/im/s, the 
KE S ~ 0.15fkj can be mainly ascribed to the swirls. On 
the other hand, when the obstacle moves with the speed 
of 150/im/s, swirls are only present at boundary of the 
SFG due to the creation of ghost vortices (vortices in low 
density regions). These ghost vortices do not contribute 
to KE S . Hence the non-zero value of KE S can be mainly 
ascribed to remnant sound waves in the superfluid. Al- 
though somewhat indistinct, the density ripples created 
by the sound are present in the central region of Fig. 

m 

In the oblate condensates, the small length of vortex 
lines along axial direction (see Fig. |4| makes them less 
susceptible to vortex bending and hence vortex reconnec- 
tions. This is due the suppression of Kelvon generation 
in oblate superfluids [48]. As a result, the vortex dipoles 
formed in pan-cake shaped traps are quite stable against 
reconnections and vortex antivortex annihilation. This 
is demonstrated in Fig. [4] here the superfluid density, 
phase, and isodensity surface with isovalue of 0.002a~^ 
are shown after 334 ms of evolution. 
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FIG. 4: (Color online) In upper row, images with labels (a) 
and (b) show the density and phase of the Fermi superfluid at 
334 ms. In lower row, half of the isosurface of the superfluid 
at 334 ms clearly shows the presence of two straight vortex 
lines. The isosurface corresponds to |0(r)| 2 = 0.002a~c. 

We wish to emphasize that the NLSE can not account 
for the presence of normal fluid in the system, and hence 
leads to the zero densities at the centers of the vortices. 
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Although very small at T = K, the normal fluid, which 
predominantly occupies the vortex cores, significantly re- 
duces the density depletion inside the vortex cores [49] - 



III. SHOCK WAVES AND SNAKE 
INSTABILITY IN BCS SUPERFLUIDS 

As mentioned in the Sec. |TJ soliton mediated decay of 
shock waves in sufficiently wide traps can also lead to the 
formation of vortex dipoles. To numerically study it, we 
create shock waves in the BCS superfluid Fermi gas by 
non-adiabatically merging the two superfluid fragments. 
The superfluid fragments are obtained by using the Gaus- 
sian obstacle potential along y-axis and z-axis, i.e., 



V ohB (r,t) = V exp{-2 { - 



x (t)Y 



(7) 



We consider Vq = 93 two, xo(t) = 0, and wq = 2.5 //m as 
the amplitude and width of the obstacle potential. The 
stationary solution with this obstacle potential is shown 
in Fig. [5|a). We achieve the non-adiabatic merging of the 
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FIG. 5: (Color online) Evolution of the superfluid density 
|0(r,t)| 2 in units of a~g C after the sudden merging of the two 
superfluid fragments. Images with labels (a), (b), (c), (d), 
(e), and (f) are the densities at t = ms, 4.8 ms, 12.7 ms, 
20.7 ms, 28.6 ms, and 36.6 ms respectively. 

two fragments of the superfluid by switching off the ob- 
stacle potential suddenly, and then letting the superfluid 
to evolve. After the obstacle potential is switched off at 
t = ms, there is the formation central peak as is shown 
in Fig. |5|b) at 4.8ms. This is consistent with the for- 
mation of central peak after the collision of two strongly 
interacting Fermi gas clouds in cigar-shaped trap [8 . The 
central peak is bounded by two band solitons which move 
outward as is shown in Fig. [5|c-d). The formation of the 
regions with large density gradients (see Fig. [5jb-c)) can 
be inferred as one of the signatures of the formation of 
shock wave. The development of large density gradients 
in superfluid, during the propagation of the shock wave, 



is also accompanied by large gradients in velocity field. 
Now, multiplying Eq. [2] by ip* and subtracting the resul- 
tant equation from its complex conjugate, one gets 



V. 



h 

Ami 







Comparing it with the continuity equation for the super- 
fluid density |^(r, t)| 2 , the velocity field of the superfluid 
is 



h (^*v^-^v^*) 



Ami 



(8) 



Using this equation, the velocity field of the superfluid 
during the initial stages of the evolution is shown in Fig. 
[6| clearly showing the development of large velocity gra- 
dients. As is evident from Fig. [6j there are larger gradi- 
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FIG. 6: (Color online) The upper and lower row show the x 
and y components of the velocity field respectively, in scaled 
units (aosc^). Images with labels (a), (b) and (c) are the x 
component of the velocity field at t — 4.8 ms, 8 ms, and 12.7 
ms respectively; exactly below them are the corresponding y 
component of the velocity field. 

ents in x component of the velocity field as compared to 
y component, indicative of the fact the potential barrier 
was along y axis. 

After some time, the soliton pair bounding the central 
peak starts to undergo snake instability [see Fig. [5jd- 
e)] leading to formation of the vortex dipoles [see Fig. 
^Jf) and[7^a)]. By this time, the outer soliton pair also 
starts showing signatures of snake instability [see Fig. 
^Jf)]. This soliton pair also decays into vortex antivortex 
pairs, Fig. [7^a-c). The three dimensional character of the 
vortex dipoles thus generated is evident from isodensity 
image in the lower panel of Fig. [7) 



IV. CONCLUSIONS 

We have studied the generation of vortex dipoles in 
weakly interacting pancake- shaped Fermi superfluid at 
T = OK using a NLSE. We find that the supercritical mo- 
tion of the repulsive obstacle across the BCS superfluid 
leads to the generation of vortex dipoles, whereas the 
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sub-critical motion can generate only sound waves. We 
also observe the decay of shock waves, which are gener- 
ated by sudden merging of the superfluid fragments, into 
vortex dipoles. The strong trapping force along axial di- 
rection provides stability to vortex dipoles against vortex 
reconnections and vortex- antivortex annihilation events. 
This is responsible for the formation the dipoles which 
can persist for several seconds in the superfluid. 
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FIG. 7: (Color online) Images with labels (a), (b), (c), are 
the superfluid densities at t = 44.6 ms, 52.5 ms, and 60.5 
ms respectively. In the lower panel, half of the isosurface of 
the superfluid at t = 52.5 ms clearly shows the presence of 
five vortex dipoles. The isosurface corresponds to |0(r)' 2 — 
0.001a" 
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